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ABSTRACT 

Following deep astrometric and photometric study of the cluster NGC 2682 (M 67), we are able to accurately determine 
its fundamental parameters. Being an old and rich cluster, M 67 is a relevant object for the analysis of the Galactic 
disk evolution. M 67 is well studied but the lack of a wide and deep Stromgren photometric study makes our results 
worthwhile. The brightest stars of the open cluster M 67 were used as uvby — Hp standard stars in our studies of 
NGC 1817 and NGC 2548, and the extension of the field covered, as well as the amount of observations, allowed to 
obtain the best set of Stromgren data ever published for this cluster. We discuss the results of our CCD uvby — Hp 
intermediate-band photometry, covering an area of about 50'x50' down to V ~ 19. Moreover, a complete membership 
segregation based on astrometric and photometric criteria is obtained. The photometric analysis of a selected sample 
of stars yields a reddening value of E(b — y) = 0.03±0.03, a distance modulus of Vo — Mv = 9.7±0.2 and [Fe/H] = 
0.01±0.14. Through isochrone fitting we found an age of logi — 9.6±0.1 (4.2±0.2 Gyr). A clump of approximately 60 
stars around V = 16, (b — y) = 0.4 could be interpreted as a population of pre-cataclysmic variable stars (if members), 
or as a stream of field G-type stars placed at twice the distance of the cluster (if non-members). 
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1. Introduction 

NGC 2682 (C0847+120), also known as M 67, in Cancer 
[a 2000 =8 h 51 m 3, 8 20 qo = +11°50'; I = 215°66, b = +31°91] 
is probably the most thoroughly studied old open cluster 
in the Galaxy, thanks to its small distance from us (esti- 
mated to be ~ 900 pc). Typically quoted values for the age 
of the cluster (~ 4 Gyr) place it among the oldest open 
clusters. In our as trometric and photometric study of sev- 
eral open clusters (jBalaguer-Nunez et al.ll2004all2005t ). the 
photometric stand ard stars were taken from this cluster 
( Nissen et al.lll987|) . thus high-quality Stromgren wide- field 
CCD photometry of it also resulted from our observations . 
The existence of a proper motions study (jZhao et alJll993f) 
of similar quality and characteristics of those used in our 
study of other clusters (based on the same plate material 
and methods, by the Shanghai Astronomical Observatory), 
allowed us to apply the same techniques and analysis tools 
to this cluster. This paper on M 67 closes the series at this 
stage of our open cluster programme. 

Photometric studies of M 67 have bee n performed by 
several authors. iMontgomerv et all (|l993f ) have presented 



a deep (V ~ 20) UBVI CCD photome try of 146 8 stars 
within 15' of the centre of the cluster, and I Fan et alj {l996) 
have studied spectrophotometry of similar depth in nine 
BATC intermediate- band filters for stars in a 1?92 x 1?92 
area centred on the cluster. In addition, variability stud- 
ies have been p erform ed in this cluster, most notably by 
iGilliland etaU (|1993|) . who conducted a very sensitive, 
highly temporally sam pled study of s t ars in the central few 
arcmin of the cluster. IStassun et al.1 (|2002f ) using differen- 
tial CCD photometry to search for variability, obtained sen- 
sitive photometry of 990 stars in a roughly square region 
one-thi rd of a degree on a side centred 5' north of the cluster 
centre. ISandauistl (|2004D also as a by-product of variability 
studies, conducted a high precision VI colour-magnitude 
diagram analysis. 



Send offprint requests to: Balaguer-Nunez, L., e-mail: 
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* Tables 3, 4 and 7 are only available in electronic form from 
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iNissen et al.1 (|l987l ) obtained accurate uvby — Hp photo- 
electric photometry of a sample of 79 stars within a radius 
of 10' from the centre, from the subgiant branch down to 
the unevolved main sequence. For this reason stars in M 67 
were used by us as photometric stand ard stars for the obser- 
vations of N GC 1817 and NGC 2548 (jBalaguer-Nunez et alj 
2004al [2005l ). The amount of observations taken in our pro- 
gramme allowed us to obtain high quality photometry of 
this cluster. In spite of being a well studied cluster, our 
Stromgren photometry deserves an analysis due both to its 
deepness and to the extension of the covered area. 
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There are al so a number of proper motion and radial 
velocity studies. [Sanders ( 1977T ) calculated probabilities of 
membership based on relative proper motions of 1866 stars 
in the cluster field with a limiting photographic magnitude 
of ~17. From 10 plate pairs of different origin and with a 
maximum ep och difference of 68 y ears, he found 649 proba- 
ble members. iGirard et al.l (|1989h gave relative proper mo- 
tions for 663 stars in an area of 42' x 34' from 44 plates with 
a maximum separation of 66 years. Although the plates 
used were of different depths, the deepest plates provided 
a limiting visual magnitude down to 16. In 1993, Zhao et 
al. derived relative proper motions for 1046 stars within a 
1?5 x 1?5 area in the region from PDS measurements of 9 
plates with a maximum epoch difference of 80 years and a 
magnitude limit V ~ 15.5. 

On the ot her hand, the radial velocity studies by 
Mathieuetal. (1986, 1990) gave precise radial- velocity mea- 
surements for 170 stars, including all main-sequence stars 
brighter than V = 12.8. 

In this paper we discuss the results of our CCD pho- 
tometric study, covering an area of about 50'x50' (Fig. |lj) 
down to V~ 19. Section 2 contains the details of the CCD 
observations and their reduction and transformation to the 
standard system. In Sect. 3 we discuss a new membership 
segregation based on the combination of parametric and 
non-parametric methods applie d to t he above-mentioned 
proper motions from Zhao et al. (|l993h . This study was se- 
lected among all other existent because it covers the biggest 
area and largest epoch difference. Also, this material is fully 
homogeneous with that used by us for the other clusters in 
this series, as from the point of view of telescope and plates, 
as from the point of view of proper motion treatment and 
reduction. This section includes a discussion on the limits 
on proper motions detection of binarity. 

In Sect. 4 we discuss the colour- magnitude diagram and 
identify the sample of probable cluster members using as- 
trometric as well as photometric criteria. Section 5 contains 
the derivation of the fundamental cluster parameters of red- 
dening, distance, metallicity and age. Section 6 studies the 
multiple star systems, blue stragglers and a new feature on 
the colour-magnitude diagram. Section 7 summarizes our 
conclusions. 



2. The Data 

2.1. Observations 

The Stromgren CCD photometry of the area was obtained 
at Calar Alto Observatory (Almeria, Spain) in January 
2000 using the 1.23 m telescope of Centro Astronomico 
Hispano-Aleman (CAHA). Further data were obtained 
at Observatorio del Roque de los Muchachos (ORM, La 
Palma, Canary Islands, Spain) in February 2000 using the 
2.5 m Isaac Newton Telescope (INT) of INC (equipped with 
the Wide-Field Camera, WFC), and in December 1998 and 
February 2000 using the 1 m Jakobus Kapteyn Telescope 
(JKT) of ING, with the Hp filters. A log of the observa- 
tions, the total number of frames, exposure times and seeing 
conditions is given in Table [1] 

We obtained photometry for a total of 1843 stars in an 
area of 50'x50' in M 67 region, down to a limiting magni- 
tude V~ 19. The area covered is shown in the finding chart 
of the cluster (Fig. [1]) . Due to the lack of Hp filter at the 
WFC-INT, it was only possible to measure it at the JKT 
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Fig. 1. Finding chart of the area under study. The cov- 
ered area is marked in black on an image of a plate 
POSSI.E -optic al R-.DSS1.486-LOW) plotted with Aladin 
Bonnarel etafll2000h . The Hp spatial coverage is limited 
to the central region. See text for details. 



and CAHA telescopes, thus limiting the spatial coverage 
with this filter. Only 288 stars (in the central region) have 
Hp values. 

2.2. Reduction and Transformation 

The reduction of the p hotome try is explained in length at 
Balaguer-Nunez et al. (|2004af ). Having the standard stars 
in the same images, the process is, in this case, equivalent 
to performe differential photometry. 

Our general procedure has been to routinely obtain twi- 
light sky flats for all the filters and a sizeable sample of 
bias frames (around 10) before and/or after every run. Flat 
fields are typically fewer in number, from five to ten per 
filter. Two or three dark frames of 2000 s were also taken. 
irafQ routines were used for the reduction process. The 
bias level was evaluated individually for each frame by av- 
eraging the counts of the most stable pixels in the overscan 
areas. The 2-D structure of the bias current was evaluated 
from the average of a number of dark frames with zero ex- 
posure time. Dark current was found to be negligible in all 
the cases. Flatfielding was performed using sigma clipped, 
median stacked, dithered twilight flats. Ten short exposures 
in every filter were taken every night with a magnitude limit 
of 19. At least one exposure per filter was pointed to 
make the centre of the cluster fall on each of the WFC 
chips. 

Our fields are not crowded. Thus, the synthetic aperture 
technique provides the most efficient measurements of rela- 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Table 1. Log of the observations used in this study. 



Telescope 


Date 


Seeing(") 


N. of frames 


Exp. Times (s) 












uvby 


Hp 


1.23 m CAHA 


2000/01/05-10 


1.1 


20 




90 


1 m JKT 


2000/02/02-06 


1.1 


42 




200 


2.5 m WFC-INT 


2000/02/02-03 


1.3 


44 


80 20 20 10 





tive fluxes within the frames and from frame to frame. We 
used the appropiate IRA F packages, and DA OPHOT and 
D AO GROW algorithms (|Stetsodll987lll990ft . We analyzed 
the magnitude growth curves and determined the aperture 
correction with the IRAF routine MKAPFILE. 

For WFC images from the INT, we employ the pipeline 
specifically developed by the Cambridge Astronomical 
Survey Unit. The process bias subtracts, gain corrects and 
flatfields the images. Catal ogues are generated using algo- 
rithms described in Irwin (1985). The pipeline gives accu- 
rate positions in right ascension and decl ination linked to 
the USN02 Catalogue (jMonet et al.lll998ft . and instrumen- 
tal magnitudes with their corresponding error s. A c omplete 
description can be found in Irwin & Lewis ( 2001) and m 
|http: //www. ast . cam. ac .uk/~wf csur /i ndex . php. 

The coefficients of the transformation equations were 
computed by a least squares method using the instrumen- 
tal magnitudes of the standard stars and the standards 
magnitudes and colours in the uvby — Hp system. Up to 
68 st andard stars in the field of the cluster (Nissen et al. 
11987ft were used depending on the size of the frame. Those 
standard stars with residuals greater than 2a were rejected. 
The reduction was performed for each night independently 
and in two steps. The first step is to determine the extinc- 
tion coefficients for each passband from the standard stars. 
With the extinction coefficients fixed, the transformation 
coefficients to the standard system were fitted. The final 
errors as a function of apparent visual magnitude are given 
in Tableland plotted in Fig. [3] 

Table 3 lists the u, v, b, y, Hp data for all 1843 stars ob- 
served in a region of 50' x 50' around the open cluster M 67 
(Fig- HI)- Star positions are given as frame {x,y) and equa- 
torial (aj2ooo Ai20oo) coordinates. An identification number 
was assigned to each star following the order of increasing 
right ascension. Column 1 is the ordinal star number; cols. 
2 and 3 are aj2ooo an d <5j20oo; cols. 4 and 5 are the respec- 
tive x, y coordinates in arcmin (arbitrary origin, close to 
cluster centre); cols. 6 and 7 are the {b — y) and its error, 
cols. 8 and 9 the V magnitude and its error, cols. 10 and 11 
the mi and its error, cols. 12 and 13 the c\ and its error, 
and cols. 14 and 15 the Hp and its error. In col. 16, stars 
considered candidate members (see Sect. 4.1.) are labelled 
'M', while those classified as non-members show the label 
'NM'. 

The cross-identification of star s in common 
with the astrometry (|Zhao et alj I1993D . BDA 
(htpp : / / ob swww . unige . ch/WEBDA) ISandersj (11977ft , 
Hipparcos (lESAl |1997|) . Tvcho-2 (|H0g et all |2000h and 
USNO-2 (jMonet et all 1 19981) catalogues is provided in 
Table 4. 




IQ 12 14 16 IS 20 
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Fig. 2. Mean internal errors of magnitude and colours (in- 
cluding the effects both from the instrumental measure- 
ments and the transformation uncertainties) as a function 
of the apparent visual magnitude, V, for all observed stars 
in the cluster region. The structure in the magnitude de- 
pendence is mainly owed to different numbers of measure- 
ments for individual stars since the centre of the cluster was 
shifted to be observed by each of the four CCD chips. 



2.3. Comparison with Previous Photometry 



After H issen et all ()l987l ). few studies of this cluster have 
been done in Stromgren photometry. I Anthony- Twaro"el 
(|1987f ) used a CCD to cover a larger are a with up to 44 star s 
in common with our sample. Comparing Ant honv-TwaroS 's 
data with ours, the mean differences in the sense ours mi- 
nus others, we get -0.01 (a = 0.04) in V, 0.00(0.02) in b~y, 
0.0 0(0.04) in mi and 0.00 (0.04) in c x . The data published 
by iJoner fc Taylor] (|1997[ ). with 14 stars in common with 
our catalogue, up to a limit of V ~ 13, give differences of 
0.01(0.03) in V, 0.01(0.0 2) in b-y, -0.02(0. 03) in mi and 
-0.01(0.04) in ci. Before INissen et aH (|l987ft . the study of 
iStrom et alj (|1971h with 19 common stars up to a limit of 
V ~ 13.5 gives: -0.02(0.03) in V, -0.03(0.02) in b - y, 
0.03(0.05) in m x , -0.02(0.06) in c x and -0.02(0.07) in Hp, 
but a large scatter and a p robable colour term in their re- 
sults was already noted by INissen et alJ (|1987ft . 

The V magnitude derived from the y filter can be com- 
pared to the published broadband data. Several studies of 
M 67 give us mean differences in V, in the sense ours minus 
others (see Table [5]). Looking at the differences, one con- 
cludes that there are no apparent systematic trends in our 
photometry. 
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Table 2. Number of stars observed (N) and mean internal errors (a) as a function of apparent visual magnitude. 



V range 




V 


(6- 


-y) 




mi 




Cl 




Hp 




N 


a 


N 


a 


N 


a 


N 


a 


N 


a 


8- 9 


3 


0.029 


3 


0.028 


3 


0.029 


2 


0.004 






9-10 


8 


0.008 


7 


0.007 


6 


0.008 


7 


0.007 


2 


0.001 


10-11 


22 


0.003 


22 


0.002 


22 


0.003 


22 


0.005 


10 


0.009 


11-12 


28 


0.004 


28 


0.022 


26 


0.003 


26 


0.005 


14 


0.016 


12-13 


114 


0.004 


114 


0.003 


114 


0.004 


114 


0.005 


56 


0.009 


13-14 


199 


0.003 


199 


0.003 


198 


0.004 


197 


0.006 


75 


0.015 


14-15 


240 


0.020 


239 


0.005 


238 


0.006 


237 


0.009 


63 


0.016 


15-16 


277 


0.007 


277 


0.010 


275 


0.011 


269 


0.018 


57 


0.019 


16-17 


299 


0.011 


299 


0.015 


272 


0.021 


233 


0.031 


10 


0.016 


17-18 


307 


0.021 


306 


0.029 


245 


0.041 


109 


0.043 


1 


0.033 


18-19 


259 


0.040 


259 


0.058 


130 


0.073 


46 


0.059 






19-20 


75 


0.093 


75 


0.125 


18 


0.141 


12 


0.106 






Total 


1831 




1828 




1547 




1274 




288 





Table 5. Comparison of the V magnitude to the published 
broadband data. 



Author 




AK>ur-oth 


a 


N 


Eegcn & Sandage (1964) 


photocl. 


-0.02 


0.04 


173 


Sanders (19891 


photocl. 


0.01 


0.05 


225 


Gilliland et al. (19931 


CCD 


-0.02 


0.03 


136 


Montgomery et al. (1993) 


CCD 


-0.01 


0.05 


967 


Kim et al. (19961 


CCD 


-0.01 


0.04 


89 


Henden (20031 


CCD 


0.01 


0.06 


306 


Sandauist (2004) 


CCD 


0.01 


0.03 


154 



3. Astrometric Analysis 

We have studied the membership segregation using para- 
metri c and non-parametric approaches i n the same way 
as in iBalaeuer-Nunez et all (l2004bl l2005h. Th anks to the 
proper motions obtained bv lZhao et al. (|l993l ) from high- 
quality plates taken with the double astrograph at the Z6- 
Se station of the Shanghai Observatory, as in the case 
of NGC 1817 and NGC 2548, we have been able to use 
the same procedure in a very homogeneous way. The fo- 
cal length of the Gautier 40 cm double astrograph at the 
Z6-Se station is 6.9 m (hence a plate scale of 30" mm -1 ). 
All the plates were measured using the Photometric Data 
Systems model 1010 automatic measuring machine at the 
Purple Mountain Observatory in Nanjing. IZhao et alj give 
relative proper motions of 1046 stars within a 1?5 x 1?5 
area, from measurements of 9 plates. The plates have a 
ma ximum epoc h diffe rence of 80 years. 

IZhao et alj (|1993t) applied the plat c-pair technique 
adopted many times at Shanghai Observatory to derive 
the relative proper motions. All linear and quadratic 
coordinate-dependent terms and the coma term were in- 
cluded in the plate-pair solution. The quoted mean errors 
of the relative proper motions vary from 0.4 mas yr _1 for 
bright stars in the inner part of the cluster field to some 
1.5 mas yr _1 for faint stars in the outer part of the clus- 
ter. T he c omparison with the proper motions of ISandeTi 
(|1977t ) and lGirard et al.1 (| 19891 ) shows a sat isfactory agree- 
ment. The magnitude limit of IZhao et all s study is V ~ 
15.5. 

There are 50 sta rs in common with the Tycho-2 

Catalogue (|Hog et al.l l2000t ) . From their absolute proper 



motions we can calculate the transformation of I Zhao et all s 
relative proper motions to the absolute reference 

frame: 

(fl a COS <5)tYC2 = 

-6.510 (±0.259) + 1.013 (±0.034) ■ (x x - 0.017 (±1.644) ■ fj, y ; 
r = 0.956; N = 44 

(M<5)tYC2 = 

-8.051 (±0.404) + 0.025 (±0.056) • [i x + 0.957 (±0.056) • 

r = 0.951; N = 50 
where r is the correlation coefficient and proper motions are 
expressed in mas yr _1 . The transformations show a rather 
good alignment of x- and y-axis with right ascension and 
declination, respectively. 

Membership determination in Zhao et al. (1993) was 
calculated with an 8-parametric Gaussian model, and a 
list of stars with probability higher than 0.8 and a dis- 
tance to the centre less than 45' gave 282 cluster members. 
For the sake of coherence with our analysis of NG C 1817 
and NGC 2548 (|Balaguer-Nunez et alJl2004bLl2005h . we ap- 
ply a 9-parametric Gaussian model and a non-paramet ric 
method to the proper motions obtained by I Zhao et al.l In 
contrast to their approach to the segregation of cluster 
members, we do not use any spatial information. Selecting 
members on kinematic and photometric criteria leaves the 
spatial information untouched, ready for a clean interpre- 
tation of radial and spatial trends. 

IZhao et al.l give x, y coordinates, relative proper mo- 
tions with their errors, number of plates used for pr oper mo- 
tion d etermination, and cross-identifications with ISandera 
(|l977l ). The authors state that their x, y coordinates are 
taken from one of the plates, but a detailed analysis of the 
data shows that this is the case for the majority of the stars, 
but there is a small subset (61 stars) whose x, y coordinates 
seem to have been taken from at least two other plates, and 
this introdu ces shifts in their quoted p ositions. Fortun ately, 
IZhao et "ail 's cross-identifications with [Sanders (1977) were 
correct, what allowed us a complete cross-identification of 
our photometric catalogue with the prop er motion d a ta not 
losing any star. This inhomogeneity in Zh ao et alj (1993) 
x, y data does not affect our results, since we do not use 
these coordinates in our study. All tests performed show 
that Zhao et al.'s proper motions are not affected. The 
mistake in the x,y coordinates of 61 stars seems to have 
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been introduced in the final preparation of the tables by 
IZhao et al.l (jl993fh 

3.1. Membership probability: general frame 

Let's consider any observational plane (a, 6), where a and 
b stand for any physical magnitudes. The distribution of 
any stellar sample on this plane can be described by its 
frequency function, 'J (a, b), measured in units of number 
of stars per unit area. The integral of the frequency func- 
tion all over the plane, i.e., its volume, gives the number of 
members of the sample, N. 

If the frequency fuction is normalised to unit vol- 
ume, then we get the probability density function (PDF), 
ip(a,b) — ^>(a,b)/N, measured in units of fractions of the 
sample per unit area. 

If two populations are superposed on the same region of 
the observational plane, but following different frequency 
functions and \I>2> then, for any point (a, b) on that 
plane, it is possible to compute the probability P for a star 
placed at that point to belong to population i (z = 1, 2): 



For our purposes, (a, b) are the proper motions and the 
observational plane is then called the vector-point diagram 
(VPD). 

3.2. The classical approach 

The parametric method assumes the existence of two pop- 
ulations in the VPD: cluster and field. The corresponding 
frequency functions are modelled as parametric Gaussian 
functions: a circular Gaussian model is adopted for the clus- 
ter distribution, while a bivariate (elliptical) Gaussian is 
admitted to describe the field. We apply a 9-pa r ametric 
Gaus sian model following iBalaguer-Nunez et alj (|2004bl 
2005) and the results are shown in Table [Sj Our implemen- 
tation takes into accou nt the individual er rors of the proper 
motion measurements (|Zhao fc Hdfl990h . If we transform 
the obtained mean relative proper motion of the cluster 
to the absolute system by the equations above, we get 
(HaCosS,^) = (-7.1±0.8,-7.6±0.4) mas yr -1 . 

3.3. The non-parametric approach 

The cluster/field segregation from astrometry has also been 
analysed with a non-parametric approach, as explained in 
length in lBalaguer-Nufiez et ail (|2004bL [2005h . We perform 
an empirical determination of the frequency functions from 
the vector-point diagram (VPD), without relying on any 
previous assumption about their profiles. In the area occu- 
pied by the cluster, the frequency function \& c +f is made 
up of two contributions: cluster, ty c , and field, \&f. To dis- 
entangle the two populations, we studied the VPD for the 
plate area outside a circle centred on the cluster. To this 
end, the centre of the cluster was chosen as the point of 
highest spatial density. The only two assumptions that we 
need to apply are (1) that it is possible to determine the 
frequency function of the field stars from some spatial area 
free of cluster members (outside the circle in our case); and, 
(2) that this frequency function is representative of the field 
frequency function in the area occupied by the cluster. 




Fig. 3. Empirical PDFs in the vector-point diagram. Top: 
ipc+f mixed sample from the inner circle of 35' centred on 
the point of maximum spatial density of stars. Centre: ipt 
field population from outside this circle. Bottom: ip c cluster 
population of M 67. It is noticeable that the cluster lies on 
a poorly populated background field. 



We did tests with circles of very different radii search- 
ing a reasonable tradeoff between cleanness (absence of 
a significant amount of cluster members) and signal-to- 
noise ratio (working area not too small, to avoid small 
number stat istics). The kernel density estimator technique 
(Hand 1982) was applied in the VPD to these data. Details 



6 



Balaguer-Nunez et al: uvby — Hp Photometry and Membership of M 67 



Table 6. Distribution parameters and their uncertainties from a 9-parametric Gaussian model applied to M 67 cluster 
and the field. n c measures the volume of the cluster frequency function (rif = 1 — n c ), fi x and fi y are the mean relative 
proper motions, er c , and a^ x and <jp_ y are the dispersions of the Gaussians fitted to the cluster and field distributions, 
and p is the correlation for the field function (it fixes the orientation of the field elliptical function). The units of /i and 
a are mas yr -1 . 





n c 




My 








P 


M 67 


0.364 


-0.59 


0.49 


0.89 










±0.016 


±0.07 


±0.06 


±0.05 








field 




0.02 


3.31 




9.00 


8.63 


-0.21 






±0.10 


±0.48 




±0.03 


±0.24 


±0.02 



of the procedure can be found in iGaladi-Enriquez et al.l 
(1998). Circular Gaussian kernel functions were used, with 
Gaussian dispersion (also called smoot hing parameter) h 
elected according to Silverman's rule (Silverman 1986). 
This way, the empirical frequency functions were computed 
for a grid with cell size of 0.2 mas yr -1 , well below the 
proper motion errors. The procedure was tested for sev- 
eral subsamples applying different proper motion cutoffs 
and the adopted one, as in the cases of NGC 1817 and 
NGC 2548, is of \fi\ < 15 mas yr- 1 . 

We finally find that the area outside a circle with a ra- 
dius of 35' centred on the cluster yields a clean frequency 
function with low cluster contamination and low noise. We 
can scale this field frequency function to represent the field 
frequency function in the inner circle, , by simply apply- 
ing a factor linked to the area. The cluster empirical fre- 
quency function can then be determined as *5f c = \I/ C _|_f — \I/f . 
Normalizing the frequency functions to volume unity we 
get the probability density functions for the mixed popu- 
lation (inside the circle), for the field (outside the circle) 
and for the cluster (non-field) population. These PDFs are 
displayed in Figure [3] We estimate the typical noise level, 
7, present in the result and we restricted the probability 
calculations to the stars with cluster PDF > 37. 

The maximum of the cluster PDF is located at 
(A*xiA*j/) = (— 0.6±0.2,0.4±0.2) mas yr -1 , which coincides 
well with the values obtained in the parametric method. If 
we transform this value to the absolute system, we get an 
absolute proper motion for the cluster of (/j a cos S, us) = 
(-7.1±0.7,-7.7±0.5) mas yr" 1 . Thirty six out of the 50 
Tycho-2 stars are classified as members giving a mean value 
of (fi a cosS,fi s ) = (-8.4±2.1,-6.1±2.2) mas yr" 1 . This is 
fully compatible with our result, but we consider our figures 
more reliable, since they are derived from the whole sample 
of stars. 



3.4. Results and discussion 

The non-parametric approach does not take into account 
the errors of the individual proper motions. But the full 
width at half maximum (FWHM) of the PDF of the clus- 
ter gives an estimation of the width of the distribution. We 
obtained a FWHM of ^4.1±0.2 mas yr -1 . Taking into ac- 
count the Gaussian dispersion owed to the smoothing pa- 
rameter h = 1.34 mas yr -1 , this would correspond to a 
r.m.s. error on pro per motions of 1.55 mas yr . But from 
Zhao et al. (l993) we know that the mean proper motion 
precision is 1.24 mas yr -1 , what gives us an intrinsic disper- 
sion component of 0.93 mas yr -1 , of the same order that 
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Fig. 4. The histogram of cluster membership probability of 
M 67. The solid line gives the results for classical parametric 
method Pp (Section l3.2p . while the dotted line corresponds 
to the non-parametric approach Pnp (Section I3.3p . The ar- 
rows mark the limiting probabilities for member selection 
for each method. 



the value obtained by the parametric membership deter- 
mination (er c = 0.89±0.05 mas yr -1 , ~ 4 km s _1 at the 
distance of 900 pc). 

Membership probabilities are computed according to 
Eq. [TJ both from the classical, parametric frequency func- 
tions (Pp) and from the non-parametric, empirical fre- 
quency funcions (Pnp). The cluster membership proba- 
bility histogram (Figure [2J shows a clear separation be- 
tween cluster members and field stars in both approaches 
(the solid line being the classical parametric method, dot- 
ted line the non-parametric approach). The non-parametric 
technique yields an expected number of cluster members 
from the integrated volume of the cluster frequency func- 
tion in the VPD areas of high cluster density, where PDF> 
37. The expected amount of cluster members in this VPD 
area is of 393. Sorting the sample in order of decreasing 
non-parametric membership probability, Pnp, the first 393 
stars are the most probable cluster members. The minimum 
value of the non-parametric probability (for the 393-rd star) 
is P NP = 0.87. 

To decide where to set the limit among members and 
non-members in the list sorted in order of decreasing para- 
metric membership probability, Pp, we accept the size of 
the cluster predicted by the non-parametric method, 393 
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stars. Thus we consider that the 393 stars of highest Pp 
are the most probable members, according to the results of 
the parametric technique. The minimum value of the para- 
metric probability (for the 393-rd star) is Pp = 0.55. It is 
worth noting that, contrary to what happened with the pre- 
vious two cluster s studied by us, and in othe r cases found 
in the literature (jGaladf-Enriquez et al.lll998f). in this c ase 
the original parametric segregation (|Zhao et al1ll993l ) at 
Pp > 0.8 was underestimating, and not overestimating, the 
number of probable members. This is because M 67 is much 
more populated, and displays a higher constrast against the 
field, than the other clusters. The need to reach a Pp cutoff 
much lower than the Pnp cutoff illustrates the effects of an 
imperfect modelling of the wings of the real distribution 
in the parametric method. Also, the importance of having 
some objective criterium to decide the Pp cutoff is evident 
from our results. 

With these limiting probabilities (Pnp > 0.87; Pp > 
0.55), we get a 97% (1014 stars) of agreement in the segre- 
gation yielded by t he two methods. Table 7 lists th e cross - 
identifications with lZhao et all (|l993l ) and lSanderi (|l977f) . 
the Pp a nd the Pnp for the 1046 stars . 

As in lBalaguer-Nunez et al.l (|2004tJ) . we accept as prob- 
able members those stars classified as member by at least 
one of the two methods. This way we get a list of 412 astro- 
metric probable member stars. Figure [5] shows the proper 
motion VPD and the sky distribution for all the measured 
stars, where empty circles denote an astrometric probable 
member of M 67, and all other stars are considered field 
stars and they are indicated by dots. 

A comparison of the cluster/field segregation for the 
143 stars in c o mmon with the radial velocity studies by 
iMathieu et all ((1986, 1990) is given in Tables [5] (candidate 
members) and [S] (non- members). The radial velocities have 
typical standard deviations for a set of measurements of 
any given star that range between 0.5 and 0.8 km s . 
To quantify the differen c es in the seg regation, we do as 
in iBalaguer-Nufiez et al.l ([2004b, 20051) and set an agree- 
ment index P c to 1 if the parametric probability, Pp , agrees 
with the radial velocity segregation, 2 if the non-parametric 
probability, Pnp, agrees, 3 if both probabilities, Pp and 
Pnp, agree and if none does. We find 125 out of 143 stars 
with P c > 0, that is 87% agreement with the radial veloc- 
ities segregation. This 13% of disagreement consists of 15 
stars out of 41 being considered non-members on the basis 
of proper motions, while only 3 out of 102 were found to 
be astrometric members but considered non-members on 
the basis of radial velocities. Those three have only one 
measurement of their radial velocities and two of them are 
known to be spectroscopic binaries (SB), therefore the mea- 
surement may be different from the systemic velocity and 
thus they cannot be completely ruled out as members. 

If we compare the parametric and non-parametric meth- 
ods, the behaviour is rather similar. For the parametric 
method we find a total of 121 stars (85%) whose member- 
ship assignation coincides with the radial velocity criterion, 
while for the non-parametric method this amounts to 124 
stars (87%). 

3.5. Magnitude effects in UCAC2 and proper motions as 
indicator of binarity 

The internal dynamics of non-resolved multiple systems 
could induce a motion of the measured photocentre dif- 



i 
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Fig. 6. UCAC2 proper motion vs magnitude for M 67 
area. Top row: ji a a nd m vs the UCAC2 magnitude 
(IZacharias et al.l 12004: botto m row: vs 2MASS J mag- 
nitude ( Skrutskie et al.l 120061 ). Proper-motion errors are 
about 1-3 mas yr _1 for stars to 12th magnitude, and 
about 4-7 mas yr _1 for fainter stars up to 16th magnitude. 
UCAC2 proper motions are taken from Tycho-2 (jHpg et al.l 
|2000| ) catalogue down to its magnitude limit V ~12. From 
this magnitude, and down to the limiting magnitude of our 
selection of candidate members (marked as empty circles) , 
there is a systematic trend depicted with a thick line 



ferent from the motion of the barycentre. This effect has 
to be null for equal-mass binaries and almost undetectable 
for binaries with components of extremely different masses, 
and should reach a maximum somewhere in between. When 
looking to the histogram of proper motion modules, this 
effect would skew the distribution inducing a longer tail 
towards high values. 

Bica fc Bonattol (|2005f ) have made an attempt of iden- 
tifying high-velocity stars as unresolved binary cluster 
members and, as a result, for mapping and quantify- 
ing the binary component in colour-magnitude diagrams. 
They an alysed 9 open cluste rs using UCAC2 proper 
motions JZacharias et al] |2004[ ) and 2MASS photometry 
(jSkrutskie et al.ll2006f) . As a test case of the method, they 
use M 67. Their study of the modulus of the proper mo- 
tion in the cluster area yields a double-peaked distribution 
instead of the skewed Gaussian that, in principle, could be 
expected. The authors associate the peaks to high-velocity 
and low-velocity populations in the cluster, and the high- 
velocity group is interpreted as produced by unresolved bi- 
nary systems. This behaviour is found in the core and in 
the halo of the cluster. The authors claim internal disper- 
sions of the order of ~6 km s" 1 for single stars and around 
~11 km s _1 for unresolved binaries, and a difference be- 
tween the two peaks of more than 20—30 km s _1 . These 
strikingly unrealistic numbers and the lack of any error es- 
timation adds suspicion to the fact that all the stars from 
the high-velocity peak are dimmer than J ~ 12 (Fig. 6 in 
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Fig. 5. The proper motion vector-point diagram (VPD) and spatial distribution of stars in M 67 region (empty circles 
for the astrometric probable members of M 67 and dots for field stars). Units of relative proper motions are mas yr _1 . 



Table 8. The cross-identi fication of stars considered as candidate members by us, with the radial velocities sudies by 
iMathieu et ail (|l986l . ll990l) and the comparison of those with the parametric (Pp) and non-parametric (Pnp) probabilities. 
See text for explanation of the agreement index P c . 



IDbda 


ID zhao 


Pp 


F 


NP 




Vr* 


IDbda 


ID zhao 


Pp 


P 


NP 


P c 


Vr" 


4 


818 





94 





96 


3 


34.0 ±2.0 


185 


958 





79 





94 


3 


33.4 ±0.0 


16 


820 





88 





96 


3 


33.2 ±0.9 


190 


989 





80 





94 





43.6 ±0.0 SB 


18 


521 





92 





95 


3 


32.9 ±0.5 


192 


603 





98 





97 


3 


32.7 ±1.1 


20 


834 





96 





97 


3 


33.1 ±1.0 


193 


999 





98 





97 


3 


33.7 ±0.5 


22 


839 





95 





96 


3 


34.4 ±1.9 SB1 


195 


1017 





97 





98 


3 


33.87±0.12 SB1 


28 


840 





97 





97 


3 


33.7 ±1.0 


210 


1005 





98 





97 


3 


33.9 ±1.1 


30 


845 





97 





97 


3 


33.6 ±0.4 


215 


608 





98 





98 


3 


33.6 ±0.6 


37 


846 





97 





97 


3 


33.7 ±0.3 


216 


1010 





97 





98 


3 


33.03±0.13 SB1 


46 


843 





98 





97 


3 


33.8 ±0.1 


217 


1002 





97 





97 


3 


33.3 ±0.4 


48 


848 





96 





98 


3 


33.0 ±0.4 


218 


1006 





98 





97 


3 


34.0 ±0.5 


51 


849 





96 





97 


3 


34.6 ±0.9 


219 


1007 





98 





98 


3 


33.4 ±0.3 SB1 


54 


863 





95 





97 


3 


34.5 ±0.7 


223 


609 





68 





94 


3 


32.8 


±0.4 


55 


862 





75 





92 





42.1 ±0.0 SB 


224 


1011 





95 





96 


3 


32.55±0.07 SB1 


72 


876 





97 





98 


3 


33.3 ±0.6 


226 


1013 





98 





98 


3 


32.5 ±1.0 


79 


860 





97 





98 


3 


34.1 ±1.7 


227 


1004 





98 





97 


3 


32.8 


±0.7 


84 


966 





78 





94 


3 


34.1 ±0.4 


231 


1029 





97 





97 


3 


32.5 


±0.6 


86 


933 





81 





86 


1 


27.6 ±0.0 


236 


1025 





98 





98 


3 


33.82±0.18 SB1 


88 


927 





96 





97 


3 


34.13±0.23 SB2 


237 


611 





97 





98 


3 


34.7 


±0.7 


95 


886 





98 





98 


3 


34.7 ±2.4 


238 


1028 





75 





94 


3 


38.3 


±0.0 


96 


929 





96 





98 


3 


33.1 ±0.7 


241 


1027 





96 





96 


3 


33.2 


±0.7 


101 


930 





97 





98 


3 


32.9 ±0.9 


243 


1026 





96 





96 


3 


33.4 ±0.5 


102 


893 





98 





97 


3 


36.17±0.17 SB1 


244 


1032 





93 





96 


3 


33.55±0.05 SB1 


104 


931 





95 





97 


3 


33.5 ±0.4 


248 


612 





98 





98 


3 


34.0 


±1.3 


105 


883 





80 





94 


3 


34.3 ±0.7 


252 


618 





97 





96 


3 


32.5 


±1.0 


108 


894 





51 





93 


2 


34.7 ±0.6 


255 


1043 





96 





97 


3 


31.3 


±0.3 


111 


891 





98 





97 


3 


33.70±0.22 SB1 


256 


226 





96 





98 


3 


33.2 


±0.9 


115 


918 





97 





98 


3 


34.4 ±0.5 


262 


1046 





97 





97 


3 


31.6 


±0.7 


117 


887 





98 





98 


3 


34.4 ±0.3 SB2 


266 


647 





78 





94 


3 


34.3 


±0.4 


119 


924 





97 





97 


3 


34.34±0.21 SB2 


271 


1058 





93 





96 


3 


34.1 


±0.8 


124 


912 





97 





98 


3 


29.6 ±0.7 SB 


272 


645 





97 





97 


3 


31.8 


±0.4 


127 


911 





97 





98 


3 


33.3 ±0.8 


281 


1060 





96 





97 


3 


33.4 


±0.8 


130 


910 





94 





98 


3 


33.7 ±0.7 


286 


697 





35 





92 


2 


33.6 


±0.5 


131 


972 





95 





96 


3 


33.3 ±2.0 SB2 


287 


698 





98 





97 


3 


32.8 


±0.5 


134 


907 





98 





97 


3 


32.5 ±1.2 SB1 


289 


683 





94 





96 


3 


32.7 


±0.4 


135 


908 





97 





97 


3 


34.3 ±0.6 


291 


700 





96 





97 


3 


34.1 


±0.1 


136 


974 





94 





96 


3 


32.87±0.12 SB1 


305 


133 





52 





92 


2 


34.2 


±0.9 


141 


947 





86 





95 


3 


33.6 ±0.4 


1182 


174 





95 





98 


3 


34.6 


±0.9 


143 


937 





91 





96 


3 


32.93±0.07 SB1 


2079 


172 





98 





98 


3 


34.4 


±0.5 


149 


949 





98 





97 


3 


35.5 ±0.8 


2087 


211 





94 





97 


3 


31.5 


±3.0 


151 


971 





81 





94 


3 


33.9 ±0.5 


3035 


1001 





98 





98 


3 


34.1 


±0.5 


157 


976 





98 





98 


3 


33.6 ±0.5 


3116 


642 





98 





97 


3 


33.60±0.13 SB1 


163 


995 





98 





98 


3 


33.8 ±0.8 


4004 


939 





98 





98 


3 


33.48±0.19 SB2 


164 


983 





79 





94 


3 


33.3 ±0.4 


4096 


578 





97 





97 


3 


34.0 ±0.8 


166 


4 





97 





98 


3 


33.3 ±0.3 


7591 


135 





91 





96 


3 


33.4 ±0.8 


170 


953 





60 





93 


3 


33.59±0.10 SB1 


7657 


505 





28 





87 


2 


33.2 ±0.9 


174 


961 





97 





98 


3 


36.2 ±4.5 


8402 


231 





83 





94 


3 


33.6 ±0.4 


176 


914 





97 





98 


3 


32.5 ±0.4 SB1 


8524 


639 





69 





89 


3 


33.6 ±0.7 


180 


985 





98 





97 


3 


35.0 ±0.8 


8571 


253 





94 





95 


3 


34.2 ±0.4 


181 


955 





98 





98 


3 


33.3 ±0.7 


8792 


282 





85 





93 


3 


33.1 ±0.5 


182 


956 





97 





98 


3 


32.2 ±2.2 


8808 


277 





85 





94 


3 


29.6 ±1.0 


184 


987 





98 





97 





61.4 ±0.0 


8832 


713 





97 





97 


3 


33.5 ±0.7 



(*) Zero errors mean that only one measurement was taken and no dispersion could be calculated. SB1: Single-lined spectroscopic binary. SB2: 
Double-lined spectroscopic binary. 



Balaguer-Nunez et al: uvby — Hp Photometry and Membership of M 67 9 
Table 9. The same as Table HI for stars considered as non-members by us. 



IDbda 


ID zhao 


Pp 


Pnp 


Pc 


Vr' 


IDbda 


ID zhao 


Pp 


Pnp 


Pc 


Vr' 


8 


822 


0.00 


0.70 


3 


29.1 ±1.8 


7116 


27 


0.00 


0.00 


3 


13.7 ±0.0 


10 


136 


0.00 


0.00 


3 


16.3 ±0.0 


7232 


52 


0.00 


0.00 


3 


4.2 ±0.0 


45 


857 


0.00 


0.00 


3 


-14.0 ±0.0 


7251 


434 


0.00 


0.47 





34.56±0.14 SB1 


49 


3 


0.00 


0.00 


3 


28.1 ±0.0 


7335 


75 


0.00 


0.00 


3 


-18.9 ±0.0 


144 


951 


0.00 


0.00 





33.72±0.14 SB1 


7350 


72 


0.00 


0.00 


3 


-13.0 ±0.0 


155 


943 


0.00 


0.00 


3 


3.6 ±1.0 


7378 


458 


0.00 


0.00 


3 


57.8 ±0.0 


173 


963 


0.00 


0.00 





33.37±0.20 SB1 


7434 


94 


0.00 


0.65 





33.7 ±0.9 


200 


601 


0.00 


0.00 


3 


12.0 ±0.4 


7440 


98 


0.00 


0.00 


3 


14.95±0.11 SB1 


206 


602 


0.00 


0.00 





34.8 ±0.3 


7445 


85 


0.00 


0.72 





33.8 ±1.0 


240 


215 


0.17 


0.84 


3 


43.3 ±0.15 SB1 


7488 


485 


0.06 


0.84 





33.0 ±0.3 


242 


1024 


0.00 


0.00 


3 


-1.1 ±0.4 


7515 


489 


0.06 


0.84 


3 


27.6 ±0.6 


277 


648 


0.00 


0.00 


3 


1.6 ±0.0 


7551 


104 


0.00 


0.67 


3 


8.9 ±0.0 


2086 


210 


0.00 


0.00 





33.3 ±0.6 


7663 


492 


0.03 


0.00 


3 


1.5 ±0.0 


3128 


620 


0.00 


0.00 


3 


-20.7 ±0.0 


7674 


496 


0.00 


0.00 





33.7 ±0.5 


4166 


481 


0.00 


0.00 


3 


-2.5 ±0.0 


8135 


566 


0.00 


0.78 





34.2 ±0.4 


4168 


482 


0.00 


0.00 


3 


39.5 ±0.4 


8355 


623 


0.00 


0.63 


3 


49.3 ±0.0 


6469 


453 


0.00 


0.00 





34.1 ±0.2 


8522 


657 


0.00 


0.00 


3 


12.6 ±0.0 


6470 


468 


0.00 


0.70 





33.3 ±0.3 


8533 


632 


0.00 


0.76 


3 


67.8 ±0.5 


6474 


495 


0.00 


0.73 


3 


9.5 ±0.6 


8557 


250 


0.00 


0.00 





33.9 ±0.5 


6514 


261 


0.00 


0.00 





34.5 ±0.7 


9015 


744 


0.10 


0.81 





33.5 ±1.6 


7112 


37 


0.00 


0.00 


3 


-9.0 ±0.0 















(*) Zero errors mean that only one measurement was taken and no dispersion could be calculated. SB1: Single-lined spectroscopic binary. SB2: 
Double-lined spectroscopic binary. 



their paper), when the abundance of binaries among the 
bright members of M 67 is well known. 

A similar study using the proper motions by I Zhao et al.l 
(1993) gives as a result a bell-shaped distribution with a tail 
towards higher velocities. We were not able to find a double- 
peaked distribution in the core nor in the halo. Different 
tries of dividing the distribution did not show any rela- 
tion to binarity when checked on the corresponding colour- 
magnitude diagram. 

As it is well depicted by the study of iPlatais et alJ 
(2003), the presence of residual colour/magnitude terms in 
any proper motion analysis is almost unavoidable. This is 
specially true when one deals with an all-sky catalogue, 
where the proper motions are derived from different first 
epoch sources. Aware of the se facts, we have chec ked for 
magnitude terms in UCAC2 (|Zacharias et al.|[2004f ) proper 
motions. We plot proper motion vs magnitude in Fig.[SJ The 
UCmag magnitude in the top row is in the UCAC bandpass 
(579-642 nm, between V and R) and should be considered 
approximate. UCAC2 proper motions are derived by us- 
ing over 140 ground- and space-based catalog ues, including 
Tych o-2 (|H0g et al.ll2000h and the AC2000.2 (|Urban et all 
119981) . In the bottom row, we have plotted the 2MASS J 
magnitude as a double check of the results. Cluster member 
candidates from our selection are marked as empty circles. 
We can note how the UCAC2 proper motions of these stars 
show a trend with magnitude for magnitudes fainter than 
from V ~ 12, highlighted by the regressi on fit marked by 
a thic k line. The behaviour detected by iBica fc Bonattol 
(|2005l ) is not related to binarity, but due to systematics 
in the proper motion data. The use of proper motion data 
from UCAC2 Catalogue should be carefully taken to avoid 
systematics. 

4. Colour-Magnitude Diagrams 

We use the V vs (v — y) colour-magnitude diagram for our 
study (Fig.[7]right and centre) because it defines the main- 
sequence of a cluster signifi cantly better t han the tradi- 
tional V vs (b — y) diagram (Meibom 2000) in presence of 
extinction. The colour-magnitude diagram of all the stars 
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Fig. 8. The colour-colour diagrams of M 67. Empty circles 
denote candidate members of M 67, chosen with astromet- 
ric and non-astrometric criteria as explained in Section 4.1. 
The thick line is the solar ZAMS standard relation shifted 
by E(b — y) = 0.03 when necessary. Evolutionary and bi- 
narity effects are noticeable in these diagrams. Some of the 
candidate members, known to be multiple stars, have very 
discordant values. 

in the area displays a fairly well defined main sequence. The 
advanced age of this cluster is obvious at first sight looking 
at the diagram. Moreover, the sequence of binaries can be 
easily followed. 

4.1. Selection of candidate member stars 

Unfortunately, proper motions are only available for the 
brightest stars in the area. Photometric measurements help 
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Fig. 7. Colour-magnitude diagrams of the M 67 area. Figure on the right shows the astrometric probable members as 
empty circles. Thick line is a shifted ZAMS, with the chosen margin for candidate members (V + 0.5,V — 1) in thin lines. 
See text for details. 



to reduce the possible field contamination in the proper 
motion membership — among bright stars — , as well as to 
enlarge the selection of members towards faint magnitudes. 
Our astrometric segregation of member stars has a limit- 
ing magnitude of V ~ 15.5. From this magnitude down to 

V ~ 18 we construct a ridge line following a fitt ing of the 
observational zero-age main sequence /ZAMS) (ICrawfordl 
[19751 H978L H971 IHilditch et al.lll983t lolsenlll984D . on the 

V vs (v — y) diagram. A selection of stars based on the dis- 
tance to this ridge line is then performed. The width of the 
sequence defined by astrometric candidate members allows 
us to establish the margins for the selection of fainter stars. 
The chosen margin includes all the stars between V + 0.5 
and V — 1 from the ridge line, as shown in the right panel 
of Fig. [7] This magnitude range allows for photometric er- 
rors in the two axis, intrinsic dispersion around the ZAMS 
and multiplicity effects. The interval of 1 magnitude at the 
bright side comprises the width of the sequence and is large 
enough to allow for equal-mass binaries (AV — 0.75) and 
for most triples (AV = 1.19 for equal-mass systems). The 
interval of 0.5 magnitudes at the faint side agrees with the 
observed width of the main sequence and is larger than the 
expected photometric errors. Luckily, the high density of 
M 67 and its outstanding contrast against the population 
of field stars (shown by the comparison of the density func- 
tions in Fig. [3]) makes the fine-tuning of these margins to be 
non-critical: the number of non-members introduced by our 
choice is low, and the loss of true members is minimised. 

Extreme outliers have been rejected making use of the 
colour-colour dia grams (Fig. Et with the help of the stan- 
dard relational dCrawfordl ll975l fl97l fl979l : IHilditch et all 
119831 ; l01senlll984f) . The astrometric and photometric selec- 
tion gives a final set of 776 candidate members. They are 
plotted in Fig.[5]as filled circles in the [mi] — [ci], m\ — (b—y), 
Ci — (b — y) and [c{\ — Hp diagrams. This photometric se- 



The asymmetry of the candidate members in the mi vsb — y 
diagram is due to: binarity, on the one hand, which in the case 
of the cool stars it always mimics a metal deficiency; and on the 
other hand, to the presence of giants, which a re located to wards 
lower mi than for dwarfs for the same b — y (|Qlsenlll984T ). 
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Fig. 9. The histograms of the distance modulus, reddening 
and metallicity of the selected member stars of M 67 with 
Hp measurements. The arrows indicate the mean values 
adopted for the cluster. 



lection is known to contain a certain amount of field con- 
tamination, reason why we prefer to talk about "candidate 
members" , rather than "bona fide members" . Anyhow, the 
selection done leads to a set of stars overwhelmingly dom- 
inated by true members: the methods applied in Sect. 5 
to derive the physical parameters of the cluster are read- 
ily capable of extracting reliable mean values from such a 
sample. 

5. Physical parameters of the cluster 

The stars of the area selected as candidate cluster members 
were classified into photometric regions and its physical pa- 
rameters dete rmine d, fol lowing the algorit hm described in 
Masana et al. (|2007f) and l Jordi et alj (|1997f ). The algorithm 
uses uvby — Hp photometry and standard relations among 
colour indices for each of the photometric regions of the HR 
diagram. 

5.1. Distance, reddening and metallicity 

Only 251 stars among the 776 candidate members have Hp 
measurements. So the computation of physical parameters 
is only possible for that subset. The results are shown in 
Fig. [9] Excluding peculiar stars and those with inconsis- 
tence among their photometric indices and applying an av- 
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erage with a 2a clipping to that subset, we found a mean 
reddening value of E(b — y) = 0.03±0.03 (corresponding 
to E(B — V) = 0.04) and a mean distance modulus of 
Vq — My = 9.4±0.4. The distance modulus may be bi- 
ased towards short values due to the presence of multi- 
ple stars treated as if they were single ones. A substantial 
fraction of binaries tend to have equal mass companions. 
These binaries would be well-separated from the cluster 
main-sequence, so they can be easily detected and removed. 
However, low mas ratio binaries would remain, and might 
systematically reduce the estimated distance since they are 
brighter and redder than single cluster members. If we set a 
conservative limit for stars with lower value than 9.4 in the 
distance modulus, we can see from the V vs b — y diagram 
that these lower values mainly correspond to stars above 
the main sequence of the cluster which are most proba- 
bly multiple stars. Excluding those we can get a value of 
Vq — My = 9.7±0.2 from the study of 66 stars that can be 
confidently considered singles. 

Metallicity is better calculated studyin g only t h e 81 F 
and G type stars in our sample following iMasanal |l994). 
We find a mean value of [Fe/H] = 0.01±0.14. 

5.2. Age 

The publication by Clem et al. (|2004D of empirically con- 
strained colour-temperature relations in the Stromgren 
system makes possible an isochrone fitting to our re- 
sults. The least unsatisfacto ry fitting is found for the 
iPietrinferni et alj (|2004|2006D trac ks. Figure HOI shows the 
isochrones from iPietrinferni et all (|2004[ ) for scaled solar 
models of solar metallicity and ages of 4.0, 4.2, 4.4 and 
4.6 Gyr for models with overshooting and of 3.4, 3.6, 
3.8 Gyr for canonical models. The addopted reddening and 
dista nce modulus are E(b — y ) = 0.03 and Vo — My = 9.7; 
from IPietrinferni et al.l (|2006l ) for a-enhanced models of so- 
lar metallicity and ages of 3.6, 3.8, 4.0 and 4.2 Gyr for mod- 
els with overshooting and of 3.0, 3.2, 3.4 Gyr for canonical 
models. For this case, the addopted reddening and distance 
modulus are E(b - y) = 0.05 and V - My = 9.7. Empty 
circles are candidate members. Canonical models give lower 
ages, a-enhanced models reproduce better the behaviour of 
the red giants clump and red giant branch, but seem to be- 
have worse in the lower giants branch and subgiants branch 
where models with overshooting seem to give a better fit. It 
seems apparent that none of the models provides a good fit 
over all the areas of the colour-magnitude diagram. We de- 
cided to adopte a compromise between these isochrones and 
give an estimation of the age of t = 4.2±0.2 Gyr (logt = 
9.62±0.02) in agreement with previous estimates by other 
outhors. 



6. Multiple Star Systems, Blue Stragglers and 
other peculiars 

6.1. Multiple Star Systems 

The high bin ary content of M 67 was already noted by 
iRacind (fl97lh : "more than half of the M 67 main-sequence 
stars b etween G2 and K5 appear to be unresolved multiple 
stars" . iMontgomerv et ahl (|1993| ) studied the distribution 
of binaries and calculated a 22% of equal-mass component 
binaries. Trying to account for the binaries with low mass 
ratios when studying the distribution of multiple stars from 
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Fig. 10. Isochrones from IPietrinferni et ahl |2004) for 
scaled solar models (top row) of solar metallicity and ages 
of 4.0, 4.2, 4.4 and 4.6 Gyr for models with overshooting 
(left panel) and ages of 3.4, 3.6, 3. 8 Gyr for canonica l mod- 
els (right panel). Isochrones from IPietrinferni et al.l ([2006) 
for alpha-enhanced models (bottom row) of solar metallic- 
ity and ages of 3.6, 3.8, 4.0 and 4.2 Gyr for models with 
overshooting (left panel) and ages of 3.0, 3.2, 3.4 Gyr for 
canonical models (right panel). 



a fiducial main sequence, they found a ratio of 38%. Due to 
the possible presence of ver y low mas s ratio binaries, this 
percentage is a lower limit. iFan et ahl (|l_996j) reconsidered 
the mass ratios and number of binaries in M 67, and using 
models leaving the mass ratio to randomly vary from zero to 
one, they found that the true binary fraction in this cluster 
depends critically on how to account for the contribution 
of low mass-ratio binaries. From the models, around a 50% 
of binaries seems a plausible scenario, with a binary mass 
ratio distribution more consistent with being random than 
double-peaked. 

A list of th e multiple star systems compiled by 
ISandauistl (|2004h and present in our photometry is shown 
in Table [TO] with our membership segregation, including 
in the l ast two rows the recently discovered eclipsing sys- 
tems bv lSandauistl (|2006D . Figure [Til shows the distribution 
of multiple star systems and blue stragglers in the colour- 
magnitude diagram. 

The stars S757, S1036 (EV Cnc), S1282 (AH Cnc) and 
ES379 (ET C nc), are WUMa contact binaries (where 
"S" refers to ISandersl's Il977l study and "ES" refers to 
lEggen fc Sandagef s ll964l study) . W UMa variables are one 
class of binary star that can produce blue stragglers after 
angular momentum loss causes the two stars to coalesce. 



6.2. Blue Stragglers 

Many studies have been done in M 67 and its rich and var- 
ied population of binaries and blue stragglers (BS). It has 
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Table 10. The cross-identification of known multiple star systems and their astrometric segregation from this w ork (last 
colum n). First column is our identification number (Table 3), second column that from BDA and third that from lSandersI 
(11977ft . 



ID 3 


IDbda 


IDs 




b-y 




V 


mi 




ci 






Hp 




M/NM 


609 


22 


821 





371±0 


001 


12.764±0.003 


0.150±0.002 





436±0 


003 








M 


629 


24 


760 





375±0 


001 


13.307±0.003 


0.180±0.002 





418±0 


004 


2 


683±0 


041 


M 


701 


55 


752 





181±0 


001 


11.310±0.003 


0.230±0.002 





841±0 


004 


2 


735±0 


006 


M 


777 


61 


757 





425±0 


004 


13.530±0.003 


0.124±0.006 





410±0 


008 


2 


640±0 


000 


M 


806 


65 


1077 





435±0 


002 


12.587±0.003 


0.173±0.003 





358±0 


004 


2 


638±0 


005 




905 


86 


1063 





630±0 


001 


13.538±0.003 


0.377±0.002 





244±0 


004 


2 


575±0 


008 


M 


910 


88 


1053 





417±0 


001 


12.245±0.003 


0.211±0.002 





375±0 


004 


2 


605±0 


008 


M 


924 


90 


975 





281±0 


001 


11.063±0.003 


0.166±0.002 





676±0 


004 


2 


666±0 


011 


NM 


962 


102 


2206 





463±0 


001 


12.396±0.003 


0.241±0.002 





346±0 


004 


2 


607±0 


014 


M 


992 


111 


986 





366±0 


001 


12.725±0.003 


0.177±0.002 





414±0 


003 


2 


624±0 


010 


M 


1005 


117 


999 





491±0 


001 


12.569±0.003 


0.254±0.002 





308±0 


003 


2 


595±0 


003 


M 


1015 


119 


1045 





384±0 


001 


12.539±0.003 


0.165±0.002 





393±0 


003 


2 


621±0 


004 


M 


1029 


123 


1070 





403±0 


004 


13.985±0.003 


0.167±0.004 





343±0 


005 


2 


601±0 


012 


M 


1033 


124 


997 





292±0 


001 


12.143±0.003 


0.172±0.002 





536±0 


004 


2 


646±0 


005 


M 


1046 


131 


1082 





266±0 


003 


11.260±0.003 


0.120±0.006 





731±0 


008 


2 


704±0 


007 


M 


1050 


134 


984 





367±0 


001 


12.259±0.003 


0.180±0.002 





409±0 


004 


2 


626±0 


000 


M 


1060 


136 


1072 





418±0 


001 


11.279±0.003 


0.160±0.002 





452±0 


004 


2 


629±0 


007 


M 


1101 


143 


1040 





545±0 


001 


11.466±0.003 


0.319±0.002 





330±0 


003 


2 


531±0 


037 


M 


1102 


144 


1000 





484±0 


001 


12.836±0.003 


0.253±0.002 





364±0 


004 


2 


596±0 


002 


NM 


1176 


161 


1036 





327±0 


004 


12.760±0.004 


0.159±0.005 





380±0 


008 


2 


660±0 


006 


M 


1206 


170 


1250 





823±0 


012 


9.769±0.009 


0.555±0.012 





372±0 


009 


2 


593±0 


009 


M 


1223 


173 


1264 





608±0 


001 


12.063±0.003 


0.428±0.002 





337±0 


004 


2 


589±0 


012 


NM 


1237 


176 


1234 





353±0 


002 


12.627±0.003 


0.165±0.021 





393±0 


046 


2 


626±0 


008 


M 


1300 


190 


1284 





166±0 


002 


10.912±0.003 


0.173±0.003 





908±0 


004 


2 


776±0 


003 


M 


1324 


195 


1242 





436±0 


002 


12.684±0.003 


0.197±0.003 





362±0 


004 


2 


632±0 


025 


M 


1352 


205 


1282 





644 ±0 


028 


13.334±0.010 


-0.086±0.030 





350±0 


028 


2 


460±0 


028 


M 


1351 


207 


1195 





264±0 


001 


12.301±0.003 


0.143±0.003 





551±0 


004 








M 


1402 


2 16 


1216 





387i0 


004 


12.6 1 3±0.004 


0. 136±0.004 





379i0 


004 


2 


728=1:0 


056 


M 


1412 


219 


1272 





385±0 


001 


12.530±0.003 


0.163±0.002 





380±0 


004 


2 


679±0 


075 


M 


1428 


224 


1221 





696±0 


001 


10.730±0.003 


0.517±0.003 





308±0 


004 


2 


713±0 


005 


M 


1488 


244 


1237 





583±0 


002 


10.741±0.003 


0.350±0.003 





350±0 


005 








M 


987 


1050 


972 





556±0 


003 


15.395±0.004 


0.277±0.004 





238±0 


007 


2 


593±0 


004 




1150 


3079 







668±0 


003 


15.788±0.003 


0.430±0.005 





139±0 


008 










1567 


3116 


1508 





364±0 


005 


12.823±0.004 


0.185±0.006 





386±0 


007 








M 


1085 


4004 


1024 





362±0 


001 


12.706±0.003 


0.170±0.002 





373±0 


003 


2 


624±0 


001 


M 


1126 


5808 


1113 





621±0 


001 


13.703±0.003 


0.321±0.003 





188±0 


006 








M 


1095 


5748 


1019 





513±0 


001 


14.338±0.003 


0.283±0.002 





241±0 


004 


2 


572±0 


015 


M 


277 


7440 


440 





800±0 


002 


8.992±0.003 


-0.062±0.003 





688±0 


004 








NM 


524 


32 


723 





535±0 


001 


15.690±0.003 


0.366±0.004 





341±0 


006 








M 


1670 


2134 


1601 





574±0 


009 


14.482±0.007 


0.290±0.011 





256±0 


015 








NM 



> 




Fig. 11. Multiple star systems (marked as crosses) and blue 
stragglers (diamonds) in a colour magnitude diagram in 
b — y (left panel) and u — y (right panel) of M 67. Empty 
circles are candidate member stars. Dots are non-member 
stars. A ZAMS is shown as reference. 



a super-blue straggler (S977), a single star with a mass 
of ~ 3 M0, more t han twice the cluster turn-off mass. 
Gillil and et alj (|l993l ) discovered two BSs with low ampli- 
tude 6 Scuti pulsations (S1280 and S1282) and evidence 
of longer perio d variations in other stars. Sandquist & 
Shetrone (|2003l hereafter SSh03) discuss t hat S968, S1066 
and S1263 could be long period variables. I Sandquist et all 



(|2003[ ) and Ivan den Berg et alj (|200lD point that S1082 is 
part of a triple system with two components being BSs, 
one in a close binary and the other one, the brightest com- 
ponent, showing some evidence of being a S Scuti star. 
The close binary is an RS CVn system and also an X- 
ray source. Other BSs detected in X-rays, S997 and S1072, 
have wide eccentric orbits whos e nature is still not under- 
stood (jvan den Berg et al.ll2004f ) . Also with an eccentric or- 
bit, S1284 is an X-ray source, too. 

We consider BS candidates those stars that, accord- 
ing to our photometry, are lo cated above and to the blu e 
of the main sequence turnoff (Stryker 1993; Bailyn 1995), 
and, at the same time, have not been classified as non- 
members in our astrometric analysis. Those member stars 
that are above the main sequence turnoff but slightly to 
the red are usually also included in the BSs list. Our selec- 
tion of BSs is then enlarged by the so-called yellow strag- 
glers (YS), astrometric member stars located between the 
turnoff and the giant branch, excluding those known to be 
unresolved binaries with fea s ible standard compo nents (e.g. 
iPortegies Zwart et al.1 Il997l ; Ivan den Beral 120011) . This al- 
lows a re-evaluation of the membership for not all, but a 
good part of the stars formerly proposed as BS members 
of M 67. This helps to clean out the list of bona fide BSs 
in this cluster. We can perform a further check on the reli- 
ability of the selection, drawing the V vs (u — y) diagram 
(see right panel of Figure [TTj) . Since red giants are faint 
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Fig. 12. Isochrones from iPietrinferni et al.l [20041 from left 
to right: and 2.5, 1.8, 1.2, 1.0, 0.7, 0.5, 0.4 Gyr for M 67. 



in UV , the photometric blends, which mimic BSs in visi- 
ble colou r-magnitude diagr ams are less problematic (see for 
example ISabbi et al.l |2004[). We find a total of 23 BS can- 
didates, listed in Table HU Two of them, S489 and S1466, 
had not been proposed as BSs before. S1466 was considered 
non-m ember by some authors (|Sanderslll977t iGirard et alj 
I1989D. but it shows X-ray emision as detected by Chandra 
(jvan den Berg et al J 120041) . 

The work by iDeng et alJ (|1999L he reafter DCLC99) , 
based in the extensive photometry by iFan et aD (|l996h . 
proposes a list of 24 BSs. Of them, 17 are astrometric mem- 
bers contained in Table El and 2 others (S977 and S2226) 
appear in our BS selection but they are not covered by our 
astrometric study. BSs S1031 and S1440 are not included 
in DCLC 99's list, but they had been also considered as BSs 
bv lAhumada fc Lapassetl (|1995[ ) and SSh03. 

The list by DCLC99 includes five stars not in Table ITT1 
The reason is that one of them (S975) has an astrometry 
incompatible with being a member, two other (S1267 and 
S1434) are astrometric members but are not covered by our 
photometry, and the other two (S145 and S277) are neither 
in our astrometric nor photometric studies. 

Other stars (S740, S821, S856 S116 5, S1183) proposed 
as BSs by lAhumada fc Lapassetl (jl995f ) are already consid- 
ered normal turnoff member stars by SSh03 and DCLC99, 
a conclusion supported by our r esults. S1947, also proposed 
bv lAhumada fc Lapassetl dl995) , is not covered by our stud- 
ies. 

The BS population of M 67 is believed to be ab- 
normally large with respe ct to other open clusters. 
lAhumada fc Lapassetl (|l995l ) give a ratio of the number of 
BSs to that of main-sequence stars within 2 magnitudes be- 
low the turnoff, -/Vbs/A^, of 30/200 for this cluster, with a 
mean of 0.025 for clusters with 6.5 < \ogt < 8.6. However, 
we find in our sample a ratio of A^bs/A^ — 23/289, similar 
to the 24/286 found by DCLC99, in closer agreement with 
the average found among other open clusters. 

We can derive photometric estimates of BS masses by 
comparing their position in the colour-magnitude diagram 
with the appropriate isochrones (see Figure fl^) from the set 



of IPietrinferni et~a l. (2004) with the corresponding metal- 
licities, reddening and distance modulus. We have plotted 
the standard isochrones from 2.5 Gyr (Mto = 1-46 Mq), 
1.8, 1.2, 1.0, 0.7, 0.5 and 0.4 Gyr (M TO = 2.70 Mq). 
The reported isochrones encompass the whole distribution 
of BSs, thus constraining the range of masses covered to 
1.5 M < M B s < 2.7 M for M 67. BSs with M BS < 
1.5 Mq, i.e. formed by the merging of low mass main- 
sequence stars, are still hidden in the main sequence. There 
are some BSs that cannot lie on the main sequence of any 
reasonable isochrone, but which can be well fit by the sub- 
giant branch sequences. These stars, mentioned above, have 
evolved to the thinning hydrogen burning shell phase and 
are rapidly moving to the base of the RGB, i.e., they are 
YS, and also some of these are candidates to be evolved blue 
stragglers (E-BS), i.e, a BS in its helium-burning phase (see 
iBellazzini et al.ll2002l and references therein). 

6.3, Anomalous clump in the colour-magnitude diagram 

The presence of a striking clump of stars around V ~ 16, 
b — y ~ 0.4, v — y ~ 1.1, can be noticed in Fig. [71 This 
accumulation of stars in the same area cannot be casual 
and deserves further consideration. We have dubbed this 
group Lola's Bunch, after one of the authors of this paper. 
There are approximately 60 stars in this clump on the V 
vs b — y diagram, where this accumulation reaches a den- 
sity more than twice that of the surrounding areas of the 
photometric space. Furthermore, the spatial distribution of 
the stars (Fig. [13]) shows a lack of them in the centre of the 
cluster. The average distance from the centre of these stars 
is r =19.94±7.85'. 

We have looked for photometry of different sources to 
check the reality of the bunch, discarding any bizarre pho- 
tometric effect in our data. We have seen that the only 
existing study covering an area wide enou gh to include the 
stars that induce this feature is the one bv lFan et all (| 19961 ) 
using the BATC system, that covers almost two square de- 
grees. 

In Fig. [T3] we can see a colour-magnitude diagram of 
the central 50'x50' area of M 67 in two different colours: 
mag3890 vs 3890-6075 and 3890-9190. Although the con- 
trast is lower than the factor two found in the Stromgren 
passbands, the bunch also appears as a significant ac- 
cumulation in the diagrams, with a local density equal 
to 1.5 times that of the surronding area in the colour- 
magnitude diagram mag3890 vs 3890-9190, and 1.75 times 
in the colour-magnitude diagram mag3890 vs 3890-6075. It 
is worth noting that the bunch is more noticeable when us- 
ing the BATC filter 3890, quite similar to Stromgren v. The 
lack of clump members at the central region of the cluster 
is also observed in BATC data. 

Stars placed in this photometric region of the colour- 
magnitude diagram are found all over the two square de- 
grees field covered by BATC data, but we have noticed 
that the bunch gives a better constrast when selecting the 
central 50'x50'. This is also approximately the area cov- 
ered by our Stromgren photometry, and coincides with the 
spatial area over which most astrometric members are con- 
centrated (Fig. [5]). 

Considering the main observational facts, we see that 
any of the two possible options (clump stars being mem- 
bers or non-members) places ourselves in front of an in- 
teresting finding. If those stars were mainly cluster mem- 
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Table 11. Our BS candidates in M 67, and thei r photometry. First column is our identification number (Table 3), 
second column that from BDA, and third that from I Sanders! (|1977t l. Last column gives the astrometric membership (void 
values for stars without astrometric membership determination). 



ID 3 


IDbda 


IDs 


b - y 


V 


mi 


ci 




Hp 


M/NM 


559 


16 


751 


0.330±0.001 


12.700±0.003 


0.172±0.002 


0.464± 


0.004 


2.914±0.029 


M 


761 


55 


752 


0.181±0.001 


11.310±0.003 


0.230±0.002 


0.841± 


0.004 


2.735±0.006 


M 


883 


81 


977 


0.344±0.004 


10.032±0.003 


0.016±0.005 


0.003± 


0.026 


2.706±0.005 




947 


95 


1005 


0.319±0.001 


12.673±0.003 


0.191±0.002 


0.450± 


0.003 


2.637±0.006 


M 


1033 


124 


997 


0.292±0.001 


12.143±0.003 


0.172±0.002 


0.536± 


0.004 


2.646±0.005 


M 


1042 


130 


2204 


0.295±0.001 


12.898±0.003 


0.168±0.002 


0.496± 


0.004 


2.671±0.004 


M 


1046 


131 


1082 


0.266±0.003 


11.260±0.003 


0.120±0.006 


0.731± 


0.008 


2.704±0.007 


M 


1050 


134 


984 


0.367±0.001 


12.259±0.003 


0.180±0.002 


0.409± 


0.004 


2.626±0.000 


M 


1060 


136 


1072 


0.418±0.001 


11.279±0.003 


0.160±0.002 


0.452± 


0.004 


2.629±0.007 


M 


1140 


153 


968 


0.064±0.001 


11.267±0.003 


0.209±0.002 


0.987± 


0.004 


2.803±0.008 


M 


1 154 


156 


1066 


n ncn_l_n nm 


in n/io_l_n nnQ 
10.94o±U.U0o 


n nnQ 
0.2U3±U.U0o 


0.997i 


0.004 


Z.84 / ±U.U0o 


M 


1176 


161 


1036 


0.327±0.004 


12.760±0.004 


0.159±0.005 


0.380± 


0.008 


2.660±0.006 


M 


1272 


184 


1280 


0.165±0.001 


12.223±0.003 


0.174±0.002 


0.820± 


0.004 


2.779±0.011 


M 


1274 


185 


1263 


0.126±0.001 


11.062±0.003 


0.218±0.002 


0.949± 


0.004 


2.639±0.012 


M 


1300 


190 


1284 


0.166±0.002 


10.912±0.003 


0.173±0.003 


0.908± 


0.004 


2.776±0.003 


M 


1351 


207 


1195 


0.264±0.001 


12.301±0.003 


0.143±0.003 


0.551± 


0.004 




M 


1370 


210 


1273 


0.375±0.001 


12.220±0.003 


0.160±0.002 


0.403± 


0.003 


2.641±0.011 


M 


1649 


282 


1440 


0.349±0.005 


13.441±0.005 


0.112±0.007 


0.442± 


0.007 




M 


1087 


4006 


1031 


0.306±0.001 


13.277±0.003 


0.171±0.002 


0.428± 


0.003 


2.648±0.008 


M 


1181 


9226 


2226 


0.294±0.004 


12.596±0.004 


0.139±0.005 


0.517± 


0.006 






1529 


261 


1466 


0.146±0.002 


10.577±0.003 


0.275±0.003 


0.796± 


0.004 




M 


261 


7489 


489 


0.321±0.002 


12.903±0.003 


0.187±0.003 


0.478± 


0.004 




M 


659 


30 


792 


0.403±0.002 


11.993±0.003 


0.132±0.003 


0.485± 


0.004 


2.732±0.014 


M 




Fig. 14. The BATC colour-colour diagrams in the central 50x50' area of M 67. The BATC filters are marked by their 
wavelenghts. 

nary systems. Second, their distribution in the outer region 
of the cluster could only be explained in terms of mass seg- 
regation, implying that the pairs have had enough time to 
experience this dynamical process since the mass-loss that 
made them as light as they would be now. 

The other option, having a clump made up from field 
stars, is no less challenging. The relatively narrow intervals 
of apparent magnitude and colour means that all these stars 
have similar spectral type and, thus, similar absolute mag- 
nitude and distance. In this case, we would have to explain 
how could be formed an accumulation of field G stars, at a 
definite distance: we would have a curtain or stream of G 
stars placed twice as far as M 67. Furthermore, their spa- 
tial distribution (that seemingly avoids the centre of M 67) 
would require some explanation. 

The scarcity of the data available make any further con- 
sideration to lay into the realm of speculation: we have al- 
ready started spectroscopic observations with the aim of 
clarifying the nature of this peculiar bunch of stars. 



bers, then their photometry could be explained if they 
were detached binaries composed by a white dwarf and 
a red dwarf, i.e: pre-cataclysmic variable stars (pre-CVs). 
Current knowledge of white dwarfs and pre-CV stars in- 
dicates that the combinations requi red to obtain the ob - 
served colours are p hysically feasible (|Bergeron et al.f l995) . 
IStassun et alJ (|2002f ) have found that five stars in this pho- 
tometric area of M 67 colour-magnitude diagram are vari- 
ables of kinds compatible with the pre-CV status, but they 
failed to perceive the relevance of this concentration of stars 
in their data. Confirming the existence of a significant pop- 
ulation of pre-CV stars in such a well studied open cluster 
(and thus, of known distance, age and metallicity) would 
have very deep implications for the understanding of the 
nature and evolution of these binary systems. First, their 
clumping in a region of the colour-magnitude diagram, and 
not on a locus following the trend that would be expected 
from the white dwarf cooling sequence, would impose tight 
restrictions to the evolutionary history of the precursor bi- 
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Fig. 13. Spatial distribution of clump stars (crosses) and 
members from our astrometric study (empty circles) in 
M 67 area. Dots represent the stars present in our CCD 
photometric sample: their spatial coverage does not reach 
the lower right nor the upper left corners. Clump stars, 
selected exclusively on a photometric basis, are disposed 
around the cluster centre, at an average distance of 20' and 
leaving a void in the middle. North is up, East is left. 

7. Conclusions 

In this paper we give a catalogue of accurate Stromgren 
uvby — Hp and J2000 coordinates for 1843 stars in an area 
of 50'x50' around M 67. We give a selection of proba- 
ble members of M 67 combining this photometric study 
with astrometric analysis using parametric as well as non- 
parametric approaches. A better determination of this clus- 
ter physical parameters based on our accurate photom- 
etry gives: E(b - y) = 0.03±0.03, [Fe/H] = 0.01±0.14, 
a distance modulus of Vb — My = 9.7±0.2 and an age 
of 4.2±0.2 Gyr (log£=9. 6±0.1). The value s are coherent 
with previous studies (see lChen et al.l (|2003f) and references 
therein) . 

An anomalous bunch of more than 60 stars in the colour- 
magnitude diagram around V — 16, (b-y) = 0.4 is tenta- 
tively interpreted as composed of pre-CVs belonging to the 
cluster, or as a stream of G stars placed twice as far as the 
cluster, two alternatives to be tested with further observa- 
tions. 



Acknowledgements. We would like to thank Simon Hodgkin and Mike 
Irwin for their inestimable help on the reduction of the images taken at 
the WFC-INT. L.B-N, also wants to thank Gerry Gilmore and Floor 
van Leeuwen for their continuous help and valuable comments, as 
well as all the people at the IoA (Cambridge) for a very pleasant stay. 
L.B-N. gratefully acknowledges financial support from EARA Marie 
Curie Training Site (EASTARGAL) during her stay at IoA. Based on 
observations made with the INT and JKT telescopes operated on the 
island of La Palma by the RGO in the Spanish Observatorio del Roque 
de Los Muchachos of the Instituto de Astrofi'sica de Canarias, and with 
the 1.52 m telescope of the Observatorio Astronomico Nacional (OAN) 
and the 1.23 m telescope at the German-Spanish Astronomical Center, 



Calar Alto, operated jointly by Max-Planck Institut fur Astronomic 
and Instituto de Astrofi'sica de Andalucia (CSIC). This study was also 
partially supported by the contracts No. AYA2003-07736, AYA2006- 
15623-C02-02, with MCYT. This research has made use of Aladin, 
developed by CDS, Strasbourg, France. This research has made use 
of the WEBDA database, operated at the Institute for Astronomy of 
the University of Vienna. 



References 

Ahumada, J. & Lapasset, E. 1995, A&AS, 109, 375 
Anthony-Twarog, B. J. 1987, AJ, 93, 647 
Bailyn, C. D. 1995, ARA&A, 33, 133 

Balaguer-Nunez, L., Jordi, C, & Galadi'-Enriquez, D. 2005, A&A, 
437, 457 

Balaguer-Nunez, L., Jordi, C, Galadi-Enriquez, D., & Masana, E. 

2004a, A&A, 426, 827 
Balaguer-Nunez, L., Jordi, C, Galadi'-Enriquez, D., & Zhao, J. L. 

2004b, A&A, 426, 819 
Bellazzini, M., Pecci, F. F., Messineo, M., Monaco, L., & Rood, R. T. 

2002, AJ, 123, 1509 
Bergeron, P., Wesemael, F., & Beauchamp, A. 1995, PASP, 107, 1047 
Bica, E. & Bonatto, C. 2005, A&A, 431, 943 

Bonnarel, F., Fernique, P., Bienayme, O., et al. 2000, A&AS, 143, 33 

Chen, L., Hou, J. L., & Wang, J. J. 2003, AJ, 125, 1397 

Clem, J. L., VandenBerg, D. A., Grundahl, F., & Bell, R. 2004, AJ, 

127, 1227 
Crawford, D. L. 1975, AJ, 80, 955 
Crawford, D. L. 1978, AJ, 83, 48 
Crawford, D. L. 1979, AJ, 84, 1858 

Deng, L., Chen, R., Liu, X., & Chen, J. S. 1999, ApJ, 524, 824 

Eggen, O. J. & Sandage, A. R. 1964, ApJ, 140, 130 

ESA. 1997, The Hipparcos and Tycho Catalogues (Noordwijk, 

Netherlands: ESA Publications Division,SP-1200) 
Fan, X., Burstein, D., & et al., J.-S. C. 1996, AJ, 112, 628 
Galadi'-Enriquez, D., Jordi, C, & Trullols, E. 1998, A&A, 337, 125 
Gilliland, R. L., Brown, T. M., McCarthy, H. K. J., et al. 1993, AJ, 

106, 2441 

Girard, T. M., Grundy, W. M., Lopez, C. E., & Altena, W. F. V. 
1989, AJ, 98, 227 

Hand, D. 1982, Kernel Discriminant Analysis (Chichester Research 
Studies Press) 

Henden, A. 2003, //ftp. nofs.navy.mil/pub/outgoing/aah/sequence/ 
Hilditch, R. W., Hill, G., & Barnes, J. V. 1983, MNRAS, 204, 241 
H0g, E., Fabricius, C, Makarov, V. V., et al. 2000, A&A, 357, 367 
Irwin, M. J. 1985, MNRAS, 214, 575 

Irwin, M. J. & Lewis, J. 2001, New Astron. Rev, 45, 105 

Joner, M. D. & Taylor, B. J. 1997, PASP, 109, 1122 

Jordi, C, Masana, E., Figueras, F., & Torra, J. 1997, A&AS, 123, 83 

Kim, S.-L., Chun, M.-Y., Park, B.-G., & Lee, S.-W. 1996, J.Korean 

Astron. Soc, 29, 43 
Masana, E. 1994, PhD thesis, Universitat de Barcelona 
Masana, E. & et al. 2007, A&A, in preparation 

Mathieu, R. D., Latham, D. W., & Griffin, R. F. 1990, AJ, 100, 1859 
Mathieu, R. D., Latham, D. W., Griffin, R. F., & Gunn, J. E. 1986, 

AJ, 92, 1100 
Meibom, S. 2000, A&A, 361, 929 

Monet, D., Bird, A., Canzian, B., et al. 1998, The USNO-A2.0 
Catalogue (U.S. Naval Observatory, Washington DC) 

Montgomery, K. A., Marschall, L. A., & Janes, K. A. 1993, AJ, 106, 
181 

Nissen, P. E., Twarog, B. A., & Crawford, D. L. 1987, AJ, 93, 634 
Olsen, E. H. 1984, A&AS, 57, 443 

Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 
168 

Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2006, ApJ, 642, 
797 

Platais, I., Wyse, R., Hebb, L., Lee, Y.-W., & Rey, S. 2003, ApJ, 591, 
L127 

Portegies Zwart, S. F., Hut, P., McMillan, S., & Verbunt, F. 1997, 

A&A, 328, 143 
Racine, R. 1971, ApJ, 168, 393 

Sabbi, E., Ferraro, F. R., Sills, A., & Rood, R. T. 2004, ApJ, 617, 
1296 

Sanders, W. L. 1977, A&AS, 27, 89 

Sanders, W. L. 1989, Revista Mex. Astron. Astrofis., 17, 31 
Sandquist, E. 2004, MNRAS, 347, 101 



16 



Balaguer-Ntmez et al: uvby — Hp Photometry and Membership of M 67 



Sandquist, E. 2006, Inf. Bull. Var. Stars, 5679 

Sandquist, E. L. & Shetrone, M. D. 2003, AJ, 125, 2173 

Sandquist, E. L., w. Latham, D., Shetrone, M. D., & Milone, A. A. E. 

2003, AJ, 125, 810 
Silverman, B. W. 1986, Density Estimation for Statistics and Data 

Analysis (J.W. Arrowsmith Ltd.) 
Skrutskic, M. F., Cutri, R. M., R. Stiening, M. W., et al. 2006, AJ, 

131, 1163 

Stassun, K. G., van den Berg, M., Mathieu, R. D., & Verbunt, F. 

2002, A&A, 382, 899 
Stetson, P. B. 1987, PASP, 99, 191 
Stetson, P. B. 1990, PASP, 102, 932 

Strom, S. E., Strom, K. M., & Bregman, J. N. 1971, PASP, 83, 768 
Stryker, L. L. 1993, PASP, 105, 1081 

Urban, S. E., Corbin, T., Wycoff, G., et al. 1998, AJ, 115, 1212 
van den Berg, M. 2001, PhD thesis, Proefschrift Universiteit Utrecht 
van den Berg, M., Orosz, J., Verbunt, F., & Stassun, K. 2001, A&A, 
375, 375 

van den Berg, M., Tagliaferri, G., Belloni, T., & Verbunt, F. 2004, 
A&A, 418, 509 

Zacharias, N., Urban, S., Zacharias, M., et al. 2004, AJ, 127, 3043 
Zhao, J. L. & He, Y. P. 1990, A&A, 237, 54 

Zhao, J. L., Tian, K. P., Pan, R. S., He, Y. P., & Shi, H. M. 1993, 
A&AS, 100, 243 



